Light cluster production in intermediate energy heavy-ion collisions induced by 

neutron-rich nuclei 
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The coalescence model based on nucleon distribution functions from an isospin-dependent trans- 
port model is used to study the production of deuteron, triton, and 3 He from heavy-ion collisions 
induced by neutron-rich nuclei at intermediate energies. It is found that the emission time of these 
light clusters depends on their masses. For clusters with same momentum per nucleon, heavier ones 
are emitted earlier. Both the yield and energy spectrum of light clusters are sensitive to the density 
dependence of nuclear symmetry energy, with more light clusters produced in the case of a stiff 
symmetry energy. On the other hand, effects due to the stiffness of the isoscalar part of nuclear 
equation of state and the medium dependence of nucleon-nucleon cross sections on light cluster 
production are unimportant. We have also studied the correlation functions of clusters, and they 
are affected by the density dependence of nuclear symmetry energy as well, with the stiff symmetry 
energy giving a stronger anti-correlation of light clusters, particularly for those with large kinetic 
energies. Dependence of light cluster production on the centrality and incident energy of heavy ion 
collisions as well as the mass of the reaction system is also investigated. 

PACS numbers: 25.70.Pq, 21.30.Fe, 21.65.+I, 24.10.Lx 



I. INTRODUCTION 



The equation of state (EOS) of an asymmetric nuclear 
matter with unequal numbers of protons and neutrons 
depends crucially on the nuclear symmetry energy. Al- 
though the nuclear symmetry energy at normal nuclear 
matter density is known to be around 30 MeV from the 
empirical liquid-drop mass formula Q, Q, its values at 
other densities are poorly known. Studies based on var- 
ious theoretical models also give widely different predic- 
tions Lack of this knowledge has hampered our un- 
derstanding of both the structure of radioactive nuclei 
S HE, | and many important issues in nuclear astro- 
physics 0, llOj , such as the nucleosynthesis during the 
pre-supernova evolution of massive stars and the prop- 
erties of neutron stars |& llPj. However, recent advances 
in radioactive nuclear beam facilities provides a unique 
opportunity to stud y th e density dependence of nuclear 
symmetry energy [3l lllllT^ . IT^ |. Theoretical studies have 
already shown that in heavy-ion collisions induced by 
neutron-rich nuclei, the effect of nuclear symmetry energy 
can be studied via the pre-equilibrium neutron/ proton 
ratio 0| , the isospin fractionation IlgL UtI fl8L , 
the isoscaling in multifragmentation 20] , the proton dif- 
ferential elliptic flow J2l| , the neutron-proton differential 
transverse flow [2^, 123} . and the tt~ to tt + ratio [24| . 
Also, it was found in a recent work that the strength of 
the correlation function for nucleon pairs with high total 
momenta in heavy-ion collisions induced by neutron-rich 
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nuclei is sensitive to the nuclear symmetry energy |25| . 
This result demonstrates that the space-time structure of 
neutron and proton emission functions is affected by the 
density dependence of nuclear symmetry energy. Since 
deuteron production in heavy-ion collisions depends on 
the space-time and momentum distributions of neutrons 
and protons at freeze out pfil ]. its production as well as 
that of other light clusters such as triton and 3 He are 
also expected to be sensitive to the density dependence 
of nuclear symmetry energy. 

To study this quantitatively, we have used in Ref. 
[27l | a coalescence model for treating cluster produc- 
tion from an isospin-dependent Boltzmann-Uchling- 
Uhlenbeck (IBUU) transport model. We found that both 
the multiplicities and energy spectra of light clusters are 
indeed sensitive to the density dependence of nuclear 
symmetry energy, with the stiff symmetry energy giv- 
ing a larger yield of clusters. In particular, the t/ 3 He 
ratio at high kinetic energies show the strongest depen- 
dence on the stiffness of nuclear symmetry energy. It was 
also pointed out in Ref. [27| that the effects due to the 
isospin-independent part of nuclear equation of state and 
the medium dependence of nucleon-nucleon cross sections 
are unimportant for light cluster production in heavy ion 
collisions induced by neutron-rich nuclei. In the present 
paper, we describe in detail the model used for study- 
ing light cluster production in heavy-ion collisions and 
present also additional results. In particular, we show 
the correlation functions of these light clusters and the 
dependence of light cluster production on the centrality 
and incident energy of heavy ion collisions as well as the 
mass of the reaction system. 

The paper is organized as follows. In Section 
ITT1 we describe the coalescence model for light clus- 
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ter production from the nucleon phase-space distribu- 
tions at freeze out. This includes the construction of 
the Wigner phase-space density functions for deuteron, 
triton, and 3 He. The isospin-dependent Boltzmann- 
Uehling-Uhlenbeck (IBUU) transport model used for 
modeling heavy ion collisions is briefly described in Sec- 
tion [ffi] Using the coalescence model given in the previ- 
ous section, the IBUU model is first applied to study light 
cluster production in heavy ion collisions involving sym- 
metric nuclei, and the results are compared with available 
experimental data. Results on light cluster production 
from heavy-ion collisions induced by neutron-rich nuclei 
at intermediate energies are presented in Section llVl We 
show not only the multiplicities and energy spectra of 
light clusters but also their dependence on the stiffness 
of nuclear symmetry energy. Dependence of light clus- 
ter production on the impact parameter, incident energy, 
and masses of the reaction system is also given. In Sec- 
tion [V] we study correlations of light clusters and their 
dependence on nuclear symmetry energy. We then con- 
clude with a summary and outlook in Section lVll 



A. The coalescence model 

In the coalescence model, the probability for produc- 
ing a cluster of nucleons is determined by the overlap of 
its Wigner phase-space density with the nucleon phase- 
space distributions at freeze out. For a system containing 
A nucleons, the momentum distribution of a M-nucleon 
cluster with Z protons can be expressed as |33| 
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where f n and f p are, respectively, the neutron and pro- 
ton phase-space distribution functions at freeze-out; p 
is the Wigner phase-space density of the M-nucleon clus- 
ter; r 4l ,-- - ,r iM _ 1 and k^,--- ,k< M _ 1 are, respectively, 
the M — 1 relative coordinates and momenta taken at 
equal time in the Af-nucleon rest frame; and G is the 
spin-isospin statistical factor for the cluster. In transport 
model simulations for heavy ion collisions, the multiplic- 
ity of a M-nucleon cluster is then given by [23] 



II. CLUSTER PRODUCTION FROM NUCLEON 
COALESCENCE 



Light cluster production has been extensively studied 
in experiments involving heavy ion collisions at all ener- 
gies (e.g., see Ref. [28j for a recent review). A popular 
model for describing the production of light clusters in 
these collisions is the coalescence model (e.g., see Ref. 
p9| for a theoretical review), which has been used at 
both intermediate [2H 0^1 |33 and high energies jU 13 . 
In most applications of the coalescence model, the en- 
ergy spectra for clusters are simply given by the prod- 
uct of the spectra of their constituent nucleons multi- 
plied by an empirical coalescence factor. In Ref. [35l| . 
for example, this simplified coalescence model based on 
the nucleon distributions obtained from the intra-nuclear 
cascade model gives a satisfactory description of mea- 
sured high energy spectra of light clusters in intermedi- 
ate energy (about 100 MeV/nucleon) heavy ion collisions. 
In more sophisticated coalescence model, the coalescence 
factor is computed from the overlap of the cluster Wigner 
phase-space density with the nucleon phase-space distri- 
butions at freeze out. The coalescence factor obtained 
from light particle emission spectra thus provides impor- 
tant information about the space-time structure of nu- 
cleon emission source and is thus a sensitive probe to the 
dynamics of heavy ion collisions [M |H S US E3, E3 . 
In the present study, this more microscopic coalescence 
model is adopted for studying light cluster production 
from heavy ion collisions induced by neutron-rich nuclei 
at intermediate energies. 
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In the above, (• • • } denotes event averaging and the sum 
runs over all possible combinations of M nucleons con- 
taining Z protons. The spin-isospin statistical factor G is 
3/8 for deuteron and 1/3 for triton or 3 He, with the lat- 
ter including the possibility of coalescence of a deuteron 
with another nucleon to form a triton or 3 He |42|. 

As pointed out in Ref. (3^ , the validity of coalescence 
model is based on the assumption that nucleon emis- 
sions are statistically independent and binding energies 
of formed clusters as well as the quantum dynamics only 
play minor roles. Since the binding energies of deuteron, 
triton and 3 He are 2.2 MeV, 7.72 MeV and 8.48 MeV, 
respectively, the coalescence model for the production of 
these light clusters in heavy-ion collisions is thus appli- 
cable if the colliding system or the emission source has 
excitation energy per nucleon or temperature above ~ 9 
MeV. Furthermore, the coalescence model is a perturba- 
tive approach and is valid only if the number of clusters 
formed in the collisions is small. As shown later, this con- 
dition is indeed satisfied for energetic deuterons, tritons, 
and 3 He in the collisions we are considering. 

In the following, we construct the Wigner phase-space 
densities for deuteron, triton, and 3 He. The method can 
be generalized, although more involved, to the Wigner 
phase-space density for four-nucleon clusters (e.g., He) 
and M-nucleon clusters. 
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B. Wigner phase-space density of deuteron 

For deuterons, we obtain their Wigner phase-space 
density from the Hulthen wave function |43j |. i.e., 



' a/3(a + (3) exp(— ar) — exp(— (3r) 
2ir(a-(3) 2 r 



(3) 



with the parameters a = 0.23 fm _1 and (3 — 1.61 fm _1 
to reproduce the measured deuteron root-mean-square 
radius of 1.96 fin |44| . The Wigner phase-space density 
of deuteron can then be obtained from 



(r+5:)^*(r-5:)expHk.R)dR, ( 4 ) 



where k = (kj^— k2)/2 is the relative momentum while 
r = (r,— Y2) is the relative coordinate of proton and neu- 
tron in the deuteron. The p^(r, k) cannot be derived 
analytically if one uses Eq. As in Ref. |3^|, we ex- 
press the Hulthen wave function in terms of the sum of 
15 Gaussian wave functions 
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where the coefficient Cj and the width parameter u>i in the 
Gaussian wave function are determined by least square 
fit, and they are given in Table D With the above 
deuteron wave function, /9^ l/ (r,k) can be obtained an- 
alytically as [HEH 
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The expansion in Eq. JSJ approximates the Hulthen 
wave function not only in coordinate space but also in 
momentum space. In momentum space, the normalized 
Hulthen wave function is given by 



4>{k) = 



y/(a + I3faf3 
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In Fig. ^ we show the exact (solid lines) and fitted 
(dashed lines) Hulthen wave function in both the coordi- 
nate and the momentum space. It is seen that the sum 
of 15 Gaussian wave functions reproduces satisfactorily 
the exact Hulthen wave function both in coordinate and 
momentum spaces. 



TABLE I: Gaussian fit coefficient d and Wi to deuteron wave 
function. 
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FIG. 1: The exact (solid curves) and fitted (dashed curves) 
Hulthen wave function in coordinate space (left panel) and 
momentum space (right panel). 



C. Wigner phase-space densities of triton and 3 He 

For triton and 3 He, their Wigner phase-space densities 
are obtained from their internal wave functions, which 
are taken to be those of a spherical harmonic oscillator 

mm, i.e., 



^(n,r 2 ,r 3 ) = (37r 2 6 4 )- 3 / 4 exp^ 



p 2 + X' 
2b 2 



In the above, we have used the normal Jacobi coordinates 
for a three-particle system, i.e., 



(9) 



where R is the center-of-mass coordinate, while p and 
A are relative coordinates. The Jacobian matrix is given 
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^From Eqs. © and Q1U[I. one obtains c?ric?r2C?r3 = 
1/| J| 3 dRdpdA =3 3 / 2 c?RdpdA. It is then easy to check 
that the wave function Eq. (|5J) is normalized to one. 
^From the relation (ri — R) + (r 2 — R) + (r 3 — R) = 
p 2 + A 2 , the root-mean-square radius r rms of triton or 
3 He can be written as 



P 2 + X 2 



\iP(r 1 ,r 2 ,v 3 )\ 2 3 3 / 2 dpdX = b 2 . (11) 



The parameter b is determined to be 1.61 fm and 1.74 
fm for triton and 3 He, respectively, from their measured 
root-mean-square radii [43] • 

The Wigner phase-space densities for triton and 3 He 
is then given by 
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where k p and k> are relative momenta, which together 
with the total momentum K arc defined by 



(13) 



with ki, k2, and kj being momenta of the three nucleons. 
The matrix J _,+ denotes the inverse of the conjugate 
complex transposition of the Jacobian matrix J, i.e., 
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their different mean-field potentials and two-body col- 
lisions in subsequent dynamic evolutions. For nucleon- 
nucleon cross sections, we use as default the experimen- 
tal values in free space. To study the effects due to the 
isospin-depcndcncc of in-medium nucleon-nucleon cross 
sections Cmedium, we also use a parameterization obtained 
from the Dirac-Brucckner approach based on the Bonn 
A potential |49j. For the experimental free-space cross 
sections, the neutron-proton cross section is about a fac- 
tor of 3 larger than the neutron-neutron or proton-proton 
cross section. On the other hand, the in-medium nucleon- 
nucleon cross sections used here have smaller values and 
weaker isospin dependence than cr eX p but strong density 
dependence. For the isoscalar potential, we use as default 
the Skyrme potential with an incompressibility Kq = 380 
MeV. This potential has been shown to reproduce the 
transverse flow data from heavy ion collisions equally well 
as a momentum-dependent soft potential with Kq = 210 
MeV|3,|E3- 

The IBUU model is solved with the test particle 
method jH^]. Although the mean-field potential is evalu- 
ated with test particles, only collisions among nucleons in 
each event are allowed. Light cluster production from co- 
alescence of nucleons is treated similar as nucleon-nucleon 
collisions, i.e., only nucleons in the same event are al- 
lowed to coalesce to light clusters. Results presented in 
the following are obtained with 20, 000 events using 100 
test particles for a physical nucleon. In the present study, 
nucleons are considered as being frozen out when their 
local densities are less than po/8 and subsequent interac- 
tions do not cause their recapture into regions of higher 
density. Other emission criteria, such as taking the nu- 
cleon emission time as its last collision time in the IBUU 
model, do not change our conclusions. All calculated re- 
sults presented in the following are in the center-of-mass 
system unless stated otherwise. 

The EOS or the energy per particle of an asymmet- 
ric nuclear matter with density p and isospin asymmetry 
6 = (p n — p p )/p, where p n and p p are, respectively, its 
neutron and proton densities, is usually approximated by 
a parabolic law [13, i-e., 
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III. IBUU MODEL AND LIGHT CLUSTER 
PRODUCTION 

To determine the space-time and momentum dis- 
tributions of nucleons at freeze out in intermediate- 
energy heavy ion collisions, we use the isospin-dependent 
Boltzmann-Uehling-Uhlenbeck (IBUU) transport model 
|H |H EH. For a review of the IBUU model, we 
refer the reader to Ref. The IBUU model includes 

explicitly the isospin degrees of freedom through differ- 
ent proton and neutron initial distributions as well as 



E(p,6)=E{p,0)+E Bym (p)6 



(15) 



where E(p, 0) is the energy per particle of symmetric nu- 
clear matter and E sym (p) is the nuclear symmetry en- 
ergy. For the density dependence of nuclear symmetry 
energy, we adopt the parameterization used in Ref. [5Jj 
for studying the properties of neutron stars, i.e., 



E sym (p) = E aym (p ) ■ u 1 , 



(16) 



where u = p/po is the reduced density and E sym (po) = 30 
MeV is the symmetry energy at normal nuclear matter 
density. In the present study, we consider the two cases 
of 7 = 0.5 (soft) and 2 (stiff) to explore the large range 
of E sym (p) predicted by many-body theories 0. 

Before investigating isospin effects on light cluster pro- 
duction from heavy-ion collisions induced by radioactive 
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FIG. 2: Experimental (open circles) and theoretical (his- 
togram) energy spectra of deuterons ((a) and (d)), tritons 
((b) and (e)), and 3 He ((c) and (f)) from 36 Ar + 58 Ni colli- 
sions at E/A = 95 MeV for impact parameter 6 = 4.5 (left 
panels) and 6.5 (right panels) fm and c .m. = 60°-120°. 



nuclei, we first compare theoretical results obtained from 
the coalescence model using nucleon phases space dis- 
tribution functions from the IBUU model with existing 
experimental data from the collisions 36 Ar + 58 Ni at 
E/A = 95 MeV at two impact parameters b = 4.5 and 6.5 
fm. In the IBUU model simulations, we use the default 
parameters, i.e., the hard EOS with Kq = 380 MeV, the 
experimental free-space cross sections <7 eX p, and the soft 
symmetry potential. However, the theoretical results are 
not corrected by experimental niters, and this may be 
partially responsible for their deviations from the exper- 
imental data. In Fig. |2 the energy spectra of deuteron, 
triton, and 3 He at center-of-mass angles c . m . = 60°-120° 
are shown in Fig. [21 by solid curves for the coalescence 
model and by open circles for the experimental data |35| . 
We see that the theoretical results reproduce satisfacto- 
rily the experimental deuteron energy spectra for both 
mid-central (b = 4.5 fm) and mid-peripheral (b = 6.5 fm) 
collisions except that they overestimate the yield at low 
kinetic energies. For the triton kinetic energy spectra, 
theoretical results also agree reasonably with experimen- 
tal data at both impact parameters except at low kinetic 
energies, where they are below the measured triton yield. 
For 3 He, both the yield and slope parameter of their 
kinetic energy spectra from theoretical calculations are 
somewhat smaller than experimental data. The flatter 
or harder energy spectra of 3 He than that of deuterons 
and tritons, known as the " 3 He puzzle" , has been seen in 
other experiments as well |H IM l5rj l57 l IH 1 159 | . 

For a more quantitative comparison, we show in Table 
ITT1 the predicted multiplicities and average kinetic ener- 



TABLE II: Calculated multiplicities and average kinetic ener- 
gies per nucleon (MeV/nucleon) of deuteron, triton, and 3 He 
together with the corresponding experimental data [35]l . 
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gies per nucleon of deuteron, triton, and 3 He in center-of- 
mass angles O c . m . = 60°-90° together with corresponding- 
experimental data. It is seen that the deuteron multiplic- 
ity from our model is about 40% above the experimen- 
tal data, while those for triton and 3 He are about 50% 
below the experimental data. As to the average kinetic 
energy per nucleon, our model describes reasonably those 
of deuteron and triton but underestimates that of 3 He. 
Failure of our model for 3 He production may be related 
to the neglect in the coalescence model of its binding en- 
ergy, which is larger than those of deuteron and triton. 
Also, the harmonic oscillator wave functions used for 3 He 
may not be as good an approximation for tritons. Fur- 
thermore, in the present dynamic coalescence model, we 
have neglected the final-state long range Coulomb po- 
tential acting on protons after they freeze out, which 
would increase the kinetic energy of produced clusters, 
particularly 3 He that has a larger charge to mass ratio 
than deuterons and tritons. In addition, the momentum- 
dependence of the isoscalar nuclear potential or symme- 
try potential, which is not considered in the present work, 
may change the particle energy spectra. Although these 
corrections need to be included for a more quantitative 
study, the symmetry energy effect on cluster production 
shown below is not expected to be appreciably affected. 



IV. ISOSPIN EFFECTS ON LIGHT CLUSTER 
PRODUCTION 

To study isospin effects on cluster production in heavy 
ion collisions, we consider mainly the reaction system of 
52 Ca + 48 Ca, which has an isospin asymmetry S = 0.2 
and can be studied at future Rare Isotope Accelerator Fa- 
cility. Another heavier reaction system of 132 Sn + 124 Sn, 
which has a similar isospin asymmetry, will also be in- 
vestigated to see how the isospin effect changes with the 
mass of the reaction system. In particular, we have con- 
sidered the collisions 52 Ca+ 48 Ca and also 132 Sn+ 124 Sn 
at an incident energy of 80 MeV/nucleon. We choose 
these reaction systems and incident energy as effects due 
to feedback from heavy fragment evaporation and feed- 
down from highly excited fragments to light cluster pro- 
duction, which cannot be described by our model, are 
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(a) Deuteron 
62 Ca+ 48 Ca 
E=80 AMeV 
b=0 fm 



expected to be less important than in collisions between 
heavier nuclei at lower energies. 



A. Yields and energy spectra of light clusters 
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FIG. 3: Time evolution of the average multiplicities of (a) 
deuterons, (b) tritons, and (c) 3 He from central collisions of 
52 Ca+ 48 Ca at E = 80 MeV/nucleon by using the soft (solid 
curves) or stiff (dashed curves) symmetry energy with a stiff 
nuclear compressibility Ko = 380 MeV and free nucleon- 
nucleon cross sections. Results using soft symmetry energy 
and free nucleon-nucleon cross sections but Ko = 201 MeV 
are shown by dotted curves, while those from soft symmetry 
energy and Ko — 380 MeV but in-medium nucleon-nucleon 
cross sections are given by dash-dotted curves. 

The yields and energy spectra of light clusters from 
heavy ion collisions are important experimental observ- 
ables that can provide useful information about particle 
production mechanism and reaction dynamics. Shown 
in Figs. |2| (a), (b) and (c) is the time evolution of the 
average multiplicities of deuteron, triton, and 3 He from 
central collisions of 52 Ca + 48 Ca at E = 80 MeV/nucleon 
by using the soft (solid curve) or stiff (dashed curve) sym- 
metry energy in the IBUU model. It is seen that produc- 
tion of these light clusters from a neutron-rich reaction 
system is sensitive to the density dependence of nuclear 
symmetry energy. The final multiplicities of deuteron, 
triton, and 3 He for the stiff symmetry energy are larger 
than those for the soft symmetry energy by 51%, 73%, 
and 100%, respectively. This is due to the fact that the 
stiff symmetry energy induces a stronger pressure in the 
reaction system and thus causes an earlier emission of 
neutrons a nd p rotons than in the case of the soft symme- 
try energy 60] , leading to a stronger correlations among 
nucleons. Furthermore, the soft symmetry energy, which 
gives a more repulsive symmetry potential for neutrons 



and more attractive one for protons in low density region 
Po) than those from the stiff symmetry energy, gen- 
erates a larger phase-space separation between neutrons 
and protons at freeze out and thus a weaker correlations 
among nucleons. The larger sensitivity of the multiplic- 
ity of 3 He to the nuclear symmetry energy than those of 
triton and deuteron seen in Fig. |3| reflects the fact that 
the effect of symmetry energy is stronger on protons than 
neutrons with lower momenta |60| . 

We have also studied the effects of using different nu- 
clear incompressibilities and nucleon-nucleon cross sec- 
tions. As shown in Fig. [21 changing the incompressibil- 
ity from K = 380 to 201 MeV (dotted curves) or using 
^medium instead of c7 GX p (dash-dotted curves) only leads 
to a small change in the yields of light clusters. This im- 
plies that the relative space-time structure of neutrons 
and protons at freeze out is not sensitive to the EOS of 
symmetric nuclear matter and the medium dependence of 
nucleon-nucleon cross sections. As explained in Ref . p5| . 
this is due to the fact that the two symmetric-matter 
EOS's create a similar pressure effect on particle emis- 
sions as they lead to different density evolutions with the 
soft one giving a larger maximum density. 
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FIG. 4: Kinetic energy spectra in center-of-mass system for 
(a) deuterons, (b) tritons, and (c) 3 He from central colli- 
sions of 52 Ca+ 48 Ca at E = 80 MeV/nucleon by using the 
soft (solid squares) or stiff (open squares) symmetry energy 
with a stiff nuclear compressibility Ko = 380 MeV and free 
nucleon-nucleon cross sections. 

The kinetic energy spectra in the center-of-mass sys- 
tem for deuterons, tritons, and 3 He's are shown in Fig. 0] 
for both the soft (solid squares) and stiff (open squares) 
symmetry energies. It is seen that the symmetry energy 
has a stronger effect on the yield of low energy clusters 
than that of high energy ones, although the effect is still 
appreciable for high energy clusters. For example, the 
symmetry energy effect on the yield of deuteron, triton, 
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and 3 He is about 60%, 100%, and 120%, respectively, if 
their kinetic energies are around 10 MeV, but is about 
30%, 40%, and 85%, respectively, if their kinetic energies 
are around 100 MeV. This follows from the fact that lower 
energy clusters are emitted later in time when the size of 
nucleon emission source is relatively independent of nu- 
clear symmetry energy, leading thus to similar probabil- 
ities for nucleons to form light clusters. Since there are 
more low energy nucleons for a stiffer symmetry energy, 
more light clusters are produced. On the other hand, 
higher energy nucleons are emitted earlier when the size 
of emission source is more sensitive to the symmetry en- 
ergy, with a smaller size for a stiffer symmetry energy. 
The probability for light cluster formation is thus larger 
for a stiffer symmetry energy. This effect is, however, re- 
duced by the smaller number of high energy nucleons if 
the symmetry energy is stiffer. As a result, production of 
high energy light clusters is less sensitive to the stiffness 
of symmetry energy. This is different from that seen in 
the correlation functions between two nucleons with low 
relative momentum, where the symmetry energy effect is 
larger for nucleon pairs with higher kinetic energies |25j| , 
as they are affected only by the size of emission source 
but not by the number of emitted nucleon pairs. 



B. The isobaric yield ratio t/ 3 He 
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FIG. 5: Time evolution of the isobaric yield ratio t/ He from 
central collisions of 52 Ca + 48 Ca at E = 80 MeV/nucleon 
using different symmetry energies, nuclear incompressibilities, 
and N-N cross sections. 

To minimize the model dependence and also other ef- 
fects, such as the feedback from heavy fragment evapora- 
tion and the feed-down from produced excited states of 
triton and 3 He, it is of interest to study the isobaric yield 
ratio t/ 3 He. In Fig. [3 we show the time evolution of the 
isobaric yield ratio t/ 3 He from central collisions of 52 Ca 
+ 48 Ca at E = 80 MeV/nucleon for different symmetry 
energies. It is seen that the stiff symmetry energy gives 
a final isobaric yield ratio t/ 3 He that is about 17% less 



than that for the soft symmetry energy. This is due to 
the fact that a soft symmetry energy increases neutron 
emission and reduces proton emission while a stiff sym- 
metry energy enhances the emission of both neutrons and 
protons. As a result, the t/ 3 He ratio is enhanced using a 
soft symmetry energy but is reduced using a stiff symme- 
try energy. It is further seen that the effect of symmetry 
energy on the t/ 3 He ratio is more pronounced at ear- 
lier stage of reactions, e.g., at 70 fm/c the t/ 3 He ratio 
is increased by 23% from the stiff to the soft symmetry 
energy. 

We have also studied the effect of nuclear incompress- 
ibility and JV-JV cross sections on the t/ 3 He ratio. As 
shown in Fig0 changing nuclear incompressibility from 
Kq = 380 to 201 MeV or using o- mc di um instead of t7 exp 
increases the ratio t/ 3 He by less than 4%. These re- 
sults imply again that the relative space-time structure 
of neutrons and protons at freeze-out is not sensitive to 
the EOS of symmetric nuclear matter and the medium 
dependence of nucleon-nucleon cross sections. 



X 




FIG. 6: The isobaric ratio t/ 3 He as a function of cluster ki- 
netic energy in the center-of-mass system by using the soft or 
stiff symmetry energy. The lines are drawn to guide the eye. 

To further investigate the symmetry energy effect, we 
show in Fig. the t/ 3 He ratio with statistical errors 
as a function of cluster kinetic energy in the center-of- 
mass system for the soft (solid squares) and stiff (open 
squares) symmetry energies. It is seen that the ratio 
t/ 3 He obtained with different symmetry energies has very 
different energy dependence. While the t/ 3 He ratio in- 
creases with kinetic energy for the soft symmetry energy, 
it decreases with kinetic energy for the stiff symmetry 
energy. The symmetry energy thus affects more strongly 
the ratio of high kinetic energy triton and 3 He. For both 
soft and stiff symmetry energies, the ratio t/ 3 He is larger 
than the neutron to proton ratio of the whole reaction 
system, i.e., N/Z= 1.5. This is in agreement with avail- 
able experimental results and also with statistical model 
simulations in other reaction systems and incident ener- 
gies nmmiiiillili!]. Itis interesting t0 note 

that although the yield of lower energy triton and 3 He 
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is more sensitive to symmetry energy than higher energy 
ones, as shown in Fig. 0] their ratio at higher energy is 
affected more by the symmetry energy. Moreover, the 
energy dependence of the ratio t/ 3 He is insensitive to 
the EOS of symmetric nuclear matter and the medium 
dependence of nucleon-nucleon cross sections. Therefore, 
the pre-equilibrium triton to 3 He ratio is a sensitive probe 
to the density dependence of nuclear symmetry energy. 



Centrality, incident energy, and reaction system 
dependence 



of smaller number of free nucleons in collisions at larger 
impact parameters. 
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FIG. 7: Impact parameter dependence of the yield for (a) 
deuterons, (b) tritons, and (c) 3 He, by using Kq — 380 MeV 
with the soft (filled squares) or stiff (open squares) symmetry 
energy. 

Light cluster production in heavy ion collisions also 
depends on the impact parameter, incident energy, and 
the mass of the reaction system. Fig. shows the im- 
pact parameter dependence of light cluster yield by us- 
ing Kq = 380 MeV with the soft (filled squares) and 
stiff (open squares) symmetry energies. It is seen that 
the symmetry energy effect on the yield of light clus- 
ters decreases slightly with increasing impact parameter. 
This simply reflects the fact that nuclear compression be- 
comes weaker when the impact parameter increases and 
the symmetry energy effect mainly comes from the dif- 
ference between the soft and stiff symmetry energies at 
low densities. Also, we find that the light cluster yield 
decreases with increasing impact parameter as a result 
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FIG. 8: Dependence of the yield for (a) deuterons, (b) tritons, 
and (c) 3 He on beam energy. The results are obtained using 
Kq = 380 MeV with the soft (filled squares) or stiff (open 
squares) symmetry energy. 

Dependence of the light cluster yield on the incident 
energy of heavy ion collisions is shown in Fig. [SJ The 
results are obtained using Kq — 380 MeV with the soft 
(filled squares) or stiff (open squares) symmetry energy. 
These results show that for the light clusters considered 
here their yields first increase and then decrease with in- 
creasing incident energy. The initial increase is due to 
increasing compression of nuclear matter with increasing 
incident energy. As a result, nucleons are emitted earlier 
and the source size for cluster production is smaller, lead- 
ing thus to a larger probability for cluster production. On 
the other hand, with further increase in incident energy, 
the phase-space volume occupied by nucleons becomes 
larger, which reduces the formation probability of light 
clusters. These features have been observed in previous 
experiments We further see that the symmetry en- 
ergy effect on light cluster production becomes weaker 
at higher incident energies. This is due to the fact that 
at higher incident energies the symmetry energy effect 
mainly comes from early compression stage of heavy ion 
collisions when the isoscalar part of nuclear equation of 
state dominates. Therefore, the colliding system expands 
rapidly at high incident energies, and the low density be- 
havior of symmetry potential thus plays a less important 
role. On the other hand, the initial compression is rela- 
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tively weak at lower incident energies, and the behavior 
of symmetry energy at low densities becomes relevant. 
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FIG. 9: Time evolution of the average multiplicity of (a) 
deuteron, (b) triton, and (c) 3 He from central collisions of 
i32 gn + i24 gn at g0 MeV/nucleon by using K = 380 MeV 
with the soft (solid curve) or stiff (dashed curve) symmetry 
energy. 

To see how the yield of light clusters changes in differ- 
ent reaction systems, we shown in Fig. |5|the time evolu- 
tion of the average multiplicities of deuteron, triton, and 
3 He from central collisions of the reaction system 132 Sn 
+ 124 Sn at 80 MeV/nucleon by using K = 380 MeV 
with the soft (solid line) or stiff (dashed line) symmetry 
energy. This reaction system has a similar isospin asym- 
metry, i.e., S — 0.22, as the reaction system 52 Ca + 48 Ca 
studied in the above, and can also be studied in future 
RIA. It is seen that the light cluster yield in 132 Sn + 
124 Sn reactions is larger than that in 52 Ca + 48 Ca reac- 
tions shown in Fig. [3] as a result of more free neutrons 
in the heavier reaction system. The symmetry energy ef- 
fect is slightly stronger in 132 Sn + 124 Sn reactions than 
in 52 Ca -I- 48 Ca. Explicitly, the final multiplicities of 
deuteron, triton, and 3 He for the stiff symmetry energy 
are larger than those for the soft symmetry energy by 
53%, 74%, and 120%, respectively. 



V. CORRELATION FUNCTIONS OF LIGHT 
CLUSTERS 

The particle emission function, which is important for 
understanding the reaction dynamics in heavy ion col- 
lisions, can be extracted from two-particle correlation 
functions; see, e.g., Refs. [H EH, IH El for reviews. 
Most extensively studied one, both experimentally and 



theoretically, has been the two-proton correlation func- 
tion |ll El El El El El Recently, data on two- 
neutron, neutron-proton correlation functions have also 
become available. The neutron-proton correlation func- 
tion is especially useful as it is free of correlations due 
to wave-function anti-symmetrization and Coulomb in- 
teractions. Indeed, Ghetti et al. have deduced from the 
measured neutron-proton correlation function the time 
sequence of neutron and proton emissions El Ell ■ The 
fragment-fragment correlation functions have also pro- 
vided useful information about the particle production 
mechanism and reaction dynamics in heayy ion collisions 
at, intermediate energies [HlIlEllIllMlIlElEl 



A. Average emission times of light clusters 
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FIG. 10: Emission time distributions for (a) protons, (b) 
deuterons, (c) tritons, and (d) 3 He by using Ko = 380 MeV 
with the soft (solid line) or stiff (dashed line) symmetry en- 
ergy. 

The emission times of different particles in heavy ion 
collisions are relevant for understanding both the colli- 
sion dynamics and the mechanism of particle production. 
In heavy ion collisions at intermediate energies, nucleon 
emissions are mainly governed by the pressure of the ex- 
cited nuclear matter during the initial stage of collisions 
[2ll |89| . Since the stiff symmetry energy gives a larger 
pressure than that due to the soft symmetry energy |60| , 
it leads to an earlier emission of neutrons and protons. To 
see the symmetry energy effect on light cluster emissions, 
we show in Figs. EK a ) ; (b), (c) and (d) the emission time 
distributions for protons, deuterons, tritons, and 3 He, re- 
spectively, by using Ko = 380 MeV with the soft (solid 
line) or stiff (dashed line) symmetry energy. One sees 
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that the emission time distributions of light clusters are 
different from that of protons. While the proton emis- 
sion time peaks earlier at ~ 45 fm/c and lasts longer, 
the emission time of light clusters peaks later at ~ 55 
fm/c and has a shorter duration. Light cluster produc- 
tion thus occurs mainly from the pre-equilibrium stage 
of heavy ion collisions as in previous BUU model simula- 
tions |83j. As to the symmetry energy effect on particle 
emissions, it is seen from Fig. ^]that the particle emis- 
sion times are earlier for the stiff symmetry energy than 
for the soft symmetry energy as the former gives a larger 
initial pressure. 
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FIG. 11: Average emission times of (a) deuterons, (b) tritons, 
and (c) 3 He as functions of their cm. kinetic energies by using 
Ko = 380 MeV with the soft (filled squares) or stiff symmetry 
energy (open squares). 

Particles emitted in earlier stage of heavy ion colli- 
sions usually have higher energy than those emitted dur- 
ing later stage of the reaction, ft is thus of interest to 
study the correlation between the emission time of a par- 
ticle and its energy. Shown in Fig. are the average 
emission times of deuterons, tritons, and 3 He as func- 
tions of their cm. kinetic energies by using K = 380 
MeV with the soft or stiff symmetry energy, ft is seen 
that the average emission time of a particle with a given 
cm. kinetic energy is earlier for the stiff symmetry en- 
ergy (open squares) than for the soft one (filled squares). 
These features are consistent with results shown in Fig. 
1101 on the particle emission time distributions. It is also 
worth noting that the emission time of a particle also 
depends on its energy, i.e., particles with higher kinetic 
energies are emitted earlier than those with lower kinetic 
energies. 

To see the effect of symmetry energy on the emission 
sequence of different particles, we show in Fig. ^| the 
average emission times of neutrons, protons, deuterons, 
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FIG. 12: Average emission times of neutrons, protons, 
deuterons, tritons, and 3 He as functions of their momenta 
per nucleon by using Ko — 380 MeV with the soft (upper 
panel) or stiff (lower panel) symmetry energy. 



tritons, and 3 He as functions of their momenta per nu- 
cleon by using Kq — 380 MeV with the soft (upper panel) 
or stiff (lower panel) symmetry energy. It is seen that or- 
dering of the average emission times for neutrons and pro- 
tons depends on the stiffness of symmetry energy, with 
a significant delay of proton emission in the case of soft 
symmetry energy. The ordering of average emission times 
for light clusters depends further on their masses, with 
heavier clusters emitted earlier than lighter ones at the 
same momentum per nucleon, which is expected from the 
coalescence model. We note that although the density de- 
pendence of nuclear symmetry energy influences signifi- 
cantly the ordering of neutron and proton average emis- 
sion times, it has only a weak effect on that for tritons and 
3 He. Experimentally, the particle emission sequence can 
be studied from the correlation functions of non-identical 
particles using various energy cuts [73, ES IH, |9(| • 



B. Correlation functions of deuteron-deuteron, 
triton-triton, and 3 He- 3 He 

In standard Koonin-Pratt formalism [9lT l92l l93l | , the 
two-particle correlation function is obtained by convolut- 
ing the emission function g(jp,x), i.e., the probability for 
emitting a particle with momentum p from the space- 
time point x = (r, t), with the relative wave function of 
the two interacting particles, i.e., 



C(P,q) 



J d 4 Xl d 4 x 2 g(P/2, ai) g (P/2, x 2 ) |<ft(q, r) 
J d 4 xi ff (P/2,xi) / d±x 2 g(P/2,x 2 ) 



(17) 
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In the above, P(= pi + P2) and q(=|(pi — P2)) are, re- 
spectively, the total and relative momenta of the particle 
pair; and 0(q, r) is the relative two-particle wave function 
with r being their relative position, i.e., r = (r 2 — ri)— 
i(vi + V2)(^2 — ti). This approach has been shown to 
be very useful in studying the reaction dynamics of inter- 
mediate energy heavy-ion collisions |68|| . In the present 
paper, we use the Koonin-Pratt method to determine the 
particle-particle correlation functions in order to study 
the effect due to the density dependence of nuclear sym- 
metry energy on the spatial and temporal structure of the 
emission source for light clusters in intermediate energy 
heavy ion collisions. 
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FIG. 13: Two-deuteron correlation functions from central col- 
lisions of 52 Ca + 48 Ca at E = 80 MeV/nucleon. Results with- 
out momentum cut are shown by solid and dashed curves for 
the soft and stiff symmetry energies, respectively. Curves with 
filled and open squares are corresponding ones for deuterons 
with total momentum P > 800 MeV/c. 

Using the program Correlation After Burner [94| . 
which takes into account final-state particle-particle in- 
teractions, we have evaluated two-cluster correlation 
functions from the emission function given by the IBUU 
model. We first show in Fig. ^| the two-deuteron cor- 
relation functions from central collisions of 52 Ca + 48 Ca 
at E — 80 MeV/nucleon. Solid and dashed curves are 
results without momentum cut using the soft and stiff 
symmetry energies, respectively. The results for deuteron 
pairs with total momenta P > 800 MeV/c are shown, 
respectively, by curves with filled and open squares for 
the soft and stiff symmetry energies. In obtaining these 
results, the final-state strong interaction between two 
deuterons includes only the relative s-wave interaction, 
which we use the one in Ref. j6^ that fits the deuteron- 
deuteron elastic scattering phase-shift by Woods-Saxon 
potential. Because of the repulsive s-wave nuclear po- 
tential and Coulomb potential, two deuterons are anti- 
correlated as shown in Fig. and the anti-correlation is 
stronger for deuteron pairs with total momenta P > 800 
MeV/c. The latter implies that deuterons with higher 
momenta are emitted earlier and thus have a smaller 



emission source size. Although the correlation function 
for all deuterons depends weakly on the stiffness of nu- 
clear symmetry energy, that for deuterons with total mo- 
menta P > 800 MeV/c clearly shows a stronger anti- 
correlation for the stiff than for the soft symmetry energy. 
For example, the soft symmetry energy gives a correla- 
tion function at q = 40 MeV/c that is 20% larger than 
that due to the stiff symmetry energy. Therefore, pre- 
equilibrium deuterons with higher momenta are emitted 
earlier or have a smaller source size for the stiff symmetry 
energy, consistent with the results shown in Fig. 1111 on 
their average emission times. The correlation function of 
deuteron pairs with high total momenta is thus sensitive 
to the density dependence of nuclear symmetry energy. 
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FIG. 14: Same as Fig. 1131 for two-triton correlation function 
with and without momentum cut. 
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FIG. 15: Same as Fig. 1131 for two- 3 He correlation function 
with and without momentum cut. 

In Figs. El an d El we show, respectively, the two- 
triton and two- 3 He correlation functions with and with- 
out momentum cut. Since the strong interaction poten- 
tial between two tritons or 3 He is poorly known, the re- 
sults obtained here include only the Coulomb potential 
in the final-state interaction as in Ref. [g^. It is seen 
that the stiff symmetry energy leads to a slightly stronger 
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anti-correlation of triton or 3 He pairs than that obtained 
with the soft symmetry energy, particularly for triton 
or 3 He pairs with total momenta P > 1000 MeV/c. The 
symmetry energy effect on heavier clusters is thus weaker 
than on lighter particles such as nucleons and deuterons, 
and this is consistent with the fact that the symmetry 
energy effect on cluster emission times becomes weaker 
with increasing particle mass as shown in Fig. 1111 

VI. SUMMARY AND OUTLOOK 

Using an isospin-dependent transport model, we have 
studied light cluster production via the coalescence model 
in heavy-ion collisions induced by neutron-rich nuclei. 
We find that the density dependence of nuclear sym- 
metry energy affects significantly the yield of deuterons, 
tritons, and 3 He in these collisions. More light clusters 
are produced for the stiff symmetry energy than the soft 
symmetry energy. This effect is particularly large when 
these clusters have lower kinetic energies. Also, the iso- 
baric ratio t/ 3 He, especially for higher energy tritons and 
3 He, shows a strong sensitivity to the density depen- 
dence of nuclear symmetry energy. It is further found 
that light cluster production is insensitive to the EOS of 
the isospin-independent part of nuclear equation of state 
and the medium dependence of nucleon-nucleon cross sec- 
tions. 

We have also studied the dependence of light cluster 
production on the impact parameter and incident energy 
of heavy ion collisions as well as the mass of the reaction 
system. We find that the symmetry energy effect on light 
cluster production depends weakly on impact parameter 
and the mass of the reaction system. However, with in- 
creasing incident energy the symmetry energy effect on 
light cluster production becomes weaker. 

We have further studied the emission time and se- 
quence as well as the correlation functions of light clus- 
ters. It is found that at same momentum per nucleon 
heavier clusters are emitted earlier than lighter ones as 
expected from the coalescence model. Since the emission 
time of light clusters is earlier for the stiff symmetry en- 
ergy than for the soft one, a stronger anti-correlation of 
light clusters is obtained with the stiff symmetry energy. 
However, the symmetry energy effect on the correlation 
functions of light clusters becomes weaker with increasing 
cluster mass. 

Isospin effects on cluster production and isotopic ratios 
in heavy ion collisions have been previously studied using 
either the lattice gas model |6^| or a hybrid of IBUU and 
statistical fragmentation model ^8|. These studies are, 
however, at lower energies than considered here, where 
effects due to multifragmcntation as a result of possible 
gas-liquid phase transition may play an important role. 



Except deuterons, both tritons and 3 He are only about 
one per event in heavy ion collisions at 80 AMeV. The 
number of other clusters such as the alpha particle is not 
large either [95l l96| . In this case, the coalescence model 
is expected to be a reasonable model for determining the 
production of light clusters from heavy ion collisions as 
demonstrated in Section IIIII This is particularly so for 
light clusters with large kinetic energies as they are not 
contaminated by contributions from decays of heavy frag- 
ments, which are mainly of low kinetic energies and may 
not be negligible in intermediate energy heavy ion colli- 
sions. Furthermore, the effect obtained in present study 
will be enhanced if other clusters are emitted earlier dur- 
ing the collisions, leading to a larger isospin asymmetry 
in the remaining emission source for light clusters. Stud- 
ies of light clusters production and correlations in heavy 
ion collisions induced by neutron-rich nuclei thus pro- 
vide another possible tool for extracting useful informa- 
tion about the density-dependence of nuclear symmetry 
energy. 

In the present work, we have not included momen- 
tum dependence in either the isoscalar mean-field poten- 
tial or the symmetry potential. The former may affect 
the properties of nuclear emission source as it has been 
shown that a momentum-dependent mean-field poten- 
tial reduces nuclear stopping or increases nuclear trans- 
parency |97| . With momentum dependence included in 
the nuclear symmetry potential, differences between the 
neutron and proton potentials are affected not only by 
the density dependence of nuclear symmetry energy but 
also by the magnitude of proton and neutron momenta 
|98l l99j | . Since the momentum-dependent symmetry po- 
tential has been shown to reduce significantly the n/p 
ratio at higher kinetic energies |100| . it is expected to 
also affect our results on the kinetic energy dependence 
of t/ 3 He ratio. It is thus of interest to study quantita- 
tively how the results obtained in the present study are 
modified by the momentum dependence of nuclear mean- 
field potentials. 
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